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TRANSIENT  WAVEFORM  ANALYSIS  OF  SWITCHING  CONVERTER 

Transient  Studies  for  a Series-Resonance 
DC  to  DC  Converter  Using  Anti-Parallel  Diodes 
and  Controlled  Transistor  Switching  Time  (Computer  Study) 


INTRODUCTION 

Advanced  military  power  processing  technology  is  guided  by  three  sig- 
nificant objectives:  1.  Reduction  in  weight  and  size  of  power  processors 

which  are  necessary  links  between  generators  of  raw  electrical  power  and 
sophisticated  user  mission  equipment;  2.  Improvement  in  efficiency  through 
reduction  in  operating  losses  which  is  an  important  key  to  thermal  design 
consideration  and  component  packaging  density,  indirectly  affecting  size 
and  weight;  3«  Component  and  operational  reliability  as  expressed  in  mean 
time  before  failure  (MTBF). 

The  present  trend  in  this  technology  clearly  indicates  emphasis  on 
advanced  switching  technology  using  semiconductor  devices  as  the  prime  power 
processors.  This  concept  permits  digitized  power  flow  control  between  input 
and  output  terminals  with  the  inherent  advantage  of  faster  control  response, 
better  waveform  resolution,  and  smaller  transformers,  inductors,  and  filters. 
All  of  these  advantages  are  directly  related  to  the  capability  of  handling 
higher  switching  frequencies.  In  this  respect,  present  limitations  are 
identified  by:  1.  The  inverse  relationship  between  power  handling  capabil- 
ity and  switching  time  of  power  transistors,  rectifiers,  and  silicon 
controlled  rectifiers  (SCR),  and  2.  The  interrelation  between  circuit  con- 
figuration of  the  power  processor  and  maximum  steady  state  and  transient 
stresses  experienced  in  its  critical  components. 

While  the  area  of  improved  semiconductor  power  switching  components  is 
of  great  concern  and  the  need  for  faster  switching  transistors  with  higher 
power  handling  capability  is  well  recognized,  it  is  the  power  switching 
circuit  concept  and  its  impact  on  technology  advancement  that  will  be 
focused  on  in  this  report. 

This  study  is  limited  to  an  analysis  of  a DC  to  DC  converter  configura- 
tion. The  most  significant  design  parameter  of  such  a converter  is  the 
ratio  of  input  and  output  voltages  and  their  operational  variation.  There 
are  two  components  of  this  ratio  which  are  of  importance  and  should  be 
distinguished:  1.  The  scaling  factor,  which  is  defined  as  the  ratio  of 

the  average  input  voltage  to  average  output  voltage,  and  2.  The  control 
ratio  range,  which  is  defined  as  maximum  input  voltage  over  minimum  output 
voltage  versus  minimum  input  voltage  over  maximum  output  voltage. 

With  the  exception  of  extreme  input/output  voltage  ratios,  scaling  is 
not  a significant  design  problem  in  the  DC  to  DC  converter  and  because  of 
the  usual  requirement  for  complete  input  to  output  isolation,  it  is  achieved 
by  a transformer  as  the  representative  circuit  element.  It  is  the  range  of 
the  control  ratio  which  defines  the  sophistication  of  a DC  to  DC  converter 
circuit.  The  concept  of  digital  power  flow  within  discrete  time  elements 
requires  that  a variable  slug  of  input  energy  be  processed  through  the 
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power  transfer  circuit  such  that  it  meets  the  average  output  power  require- 
ments. From  this  discussion,  it  is  obvious  that  the  range  of  change  in 
discrete  processing  time  is  directly  proportional  to  the  previously  defined 
control  ratio  range.  In  turn,  the  need  for  storage  capability  in  the  filter 
sections  (input  and  output)  is  also  proportional  to  this  range  of  discrete 
processing  time  which  is  directly  reflected  in  the  size  and  weight  of  power 
processing  circuit  components  and  filters.  At  a first  glance,  the  nature 
of  DC  input  and  output  requirements  would  make  it  ideal  to  process  discrete 
time  energy  slugs  with  rectangular  pulse  current  and  voltage  wave  shapes 
(pulse  width  modulation) . However,  achievable  dv/dt  and  di/dt  parameters 
directly  depend  on  switching  component  capability  and  circuit  design. 
Therefore,  with  the  desired  increase  in  switching  frequency,  it  is  realized 
that  the  concessions  that  have  to  be  made  in  deviating  from  ideal  rectangu- 
lar pulse  waveforms  become  so  severe  that  this  idealization  is  no  longer  a 
useful  concept.  Consequently,  a new  concept  approach  follows  logically. 

If  idealized  rectangular  wave  shapes  can  not  be  processed,  then  it  may  be 
feasible  to  go  to  sinusoidal  wave  shapes  which  are  ingredients  of  idealized 
resonance  circuits.  This  concept  has  been  widely  explored  on  an  empirical 
basis  in  the  past  five  years  and  with  remarkable  success.  One  of  the  most 
salient  features  of  this  approach  is  the  use  of  zero  current  switching 
without  the  need  of  forced  commutation,  which  considerably  simplifies  and 
reduces  circuit  requirements  and  increases  the  power  handling  rating  of  the 
switching  components.  This  in  turn  allows  one  tq  go  to  higher  switching 
frequencies  which  more  than  compensates  for  the  increased  rating  of  energy 
storage  components  that  is  related  to  the  difference  between  sinusoidal  and 
rectangular  wave  shapes. 

The  series  resonance  converter  circuitry  selected  for  analysis  in  this 
report  represents  one  of  the  simplest  circuit  configurations  of  a DC  to  DC 
switching  converter  with  a minimum  of  components  in  the  power  processing 
circuit.  It  is  true,  that  the  deviation  from  the  ideal  digitized  energy 
flow  concept  with  rectangular  wave  shapes  presents  a more  complex  problem 
to  the  analysis  of  power  transfer  in  this"  circuit.  This  has  generated  the 
need  for  a more  detailed  transient  waveform  analysis  and  a better  under- 
standing of  the  control  concept  tnat  relates  discrete  time  control  to  the 
significant  output  performance  parameters  in  the  resonance  type  converter. 1 

ANALYSIS 

The  analysis  applies  to  the  series  resonance  DC  to  DC  converter  cir- 
cuit, as  shown  in  Figure  1,  which  has  been  treated  in  some  detail  in  the 
literature. 2 The  input  symmetry  of  the  circuit,  with  respect  to  the  split 
resonance  capacitors  Cl  and  C2  and  the  two  transistor  switching  elements 
TR1  and  TR2,  permits  simplifications  in  the  equivalent  power  transfer  circuits 
which  are  applicable  for  the  4-cycle  consecutive  switching  states  as  shown 
in  Figure  2.  It  is  seen  that  these  circuits  all  contain  the  same  elements 


F.  C.  Schwarz,  "Engineering  Information  on  an  Analog  Signal  to 
Discrete  Time  Interval  Converter  (ASDTIC),"  NASA  CR-134544  Report, 
Contract  NAS  3-16 791  (NASA),  Power  Electronics  Associates,  June  1974. 

2p.  C.  Schwarz,  "An  Improved  Method  of  Resonant  Current  Pulse  Modula- 
tion for  Power  Converters,"  Proc.  IEEE  Power  Electronics  Specialists 
Conference,  PESC  '75  Record,  pp.  194-204  (1975)- 
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if  Cl  = C2;  TR1  and  TR2  are  considered  as  equivalent,  idealized  switches. 

The  transformer  coupling  of  the  output  loop  containing  a single  filter 
capacitor  (see  Figure  1)  is  represented  in  Figure  2 by  the  impressed  voltage 
Vip  on  the  primary  winding  of  this  ideal  transformer.  The  switching  states 
are  characterized  by  their  respective  sets  of  boundary  conditions  at  the 
beginning  of  each  switching  period.  These  boundary  conditions  directly 
follow  from  the  end  of  period  state  of  the  preceding  switching  cycle. 

Therefore,  the  solution  in  the  time  domain  for  any  of  the  selected  voltages 
and  currents,  which  have  been  considered  in  the  state  equations  of  the 
circuit,  represent  a continuous  wave  train  which  is  an  exact  replica  of 
waveforms  that  may  be  examined  experimentally  on  a scope,  using  a converter 
circuit  for  which  Figure  1 is  a representative  model. 

A computer  program  has  been  written  that  calculates  and  plots  out,  on 
a X-Y  recorder,  the  time  domain  solutions  of  each  state  variable  as  a 
continuous  string  of  waveforms.  This  allows  evaluation  of  transient 
stresses  and  steady  state  conditions.  Step  function  changes  of  operational 
conditions  can  be  introduced  as  a new  set  of  input  parameters  which  affect 
the  boundary  conditions  at  each  beginning  of  a switching  state.  This  allows 
studying  of  transient  and  stability  problems  and  even  permits  control  loop 
design  evaluations  on  the  basis  of  Analog  Signal  to  Discrete  Time  Interval 
Converter  (ASHTIC)  type  control  functions. 3 

The  termination  of  each  transistor  conduction  period  is  based  on  a 
test  criterion  of  zero  current  flow.  Then  the  analysis  of  the  next  diode 
conduction  period  is  executed  by  the  computer  for  the  predetermined  diode 
conduction  time  which  is  introduced  as  a dimensionless  parameter  relating 
diode  conduction  time  to  the  half -period  time  at  resonance  frequency. 

Calculated  output  parameters  are  also  dimensionless  and  normalized  with 
respect  to  the  input  voltage.  The  chosen  range  of  the  normalized  diode 
conduction  time  is  zero  to  unity  (0<Z2cl).  This  avoids  interruptions 
of  power  flow  in  the  input  circuit  for  normal  load  conditions.  Results 
of  the  analysis  are  presented  for  a converter  circuit  configuratioo..¥hich 
operates  over  a frequency  range  from  13. 5 to  27  kHz.  For  the  specified 
range  of  diode  conduction  period  of  (0<Z2«1),  a range  of  the  output  to 
input  voltage  control  ratio  (Vq/E)  of  2.5  to  0.5  can  be  utilized  at 
nominal  load. 

ASSUMPTIONS  AND  CONSTRAINTS 

j 

The  simplifying  assumptions  which  have  been  selected  for  the  analysis  .1 

are:  .11 

1.  All  passive  circuit  components  are  represented  by  lumped  } 

parameters  of  L,  C,  and  R. 

2.  Inductors  and  capacitors  are  considered  ideal  and  lossless. 


F.  C.  Schwarz,  NASA  CH-l^h^hh  Report,  op.  cit.,  p.  2. 
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3.  The  coupling  transformer  is  considered  an  ideal  voltage 
scaling  device,  without  losses  and  represented  by  a turns  ratio  M = W1/W2. 

4.  Transistors  and  diodes  are  considered  lossless,  instantaneous 
switches  characterized  by  the  two  states  on  and  off.  They  are  unipolar 
conducting  devices  requiring  the  ending  of  a discrete  switching  period 
whenever  the  current  reaches  zero.  This  requirement  determines  the  four 
equivalent  circuit  configurations  shown  in  Figure  2. 

5.  The  switching  effect  of  the  rectifier  bridge  in  the  DC  output 
circuit  is  modeled  by  a function  VT  = Sign  (lL)  • M * VQ,  implied  by  the 
change  of  sign  of  VT  whenever  the  current  changes  direction. 

6.  Continuity  of  capacitor  voltages  and  inductor  currents  is  a 
requirement  for  time  boundary  conditions  between  switching  states. 

STATE-VARIABLES  

In  each  of  the  four  switching  states,  the  circuit  state  can  be  de- 
scribed by  three  selected  state  variables.  To  facilitate  a dimensionless 
analysis,  the  state  variables  selected  are  Vq/E,  E,  and  Vp/E.  Two  other 
variables  important  for  circuit  evaluation  are  J/E,  the  normalized  value  of 
the  inductor  current,  and  Vq/E  the  normalized  value  of  the  DC  output  volt- 
age. Both  are  derived  by  simple  relations  with  Vc  and  Vp,  respectively, 
and  were  not  treated  as  independent  ■variables,  in  order  to  keep  the  coeffi- 
cient matrix  to  a simple  manageable  3 Ly  4 matrix. 

General  Expressions  for  the  State  Equations  in  S-Space 


vy,  x (3) 


A,  A ^ A • S + A , 
yxl  + yx2  + yx3  yx4 

So  . s-sxl  + sx32 


(1) 


l 


where 

y subscript  defines  the  state  variable  (Y  = 1 to  3) 

x subscript  defines  the  order  of  the  system  (x  = 0 to  3) 


The  characteristic  equation  in  S-Space  is: 


+ B^S3  + 


B^S  = 0 


(2) 


The  roots  of  Equation  (2)  are  only  a function  of  the  circuit  element 
parameters  C,  L,  RL,  Cq,  and  M.  The  values  A^  in  Equation  (l)  depend  on 
the  initial  boundary  conditions  for  each  switching  state  and  are  expressed 
by  matrix  equations  of  the  form. 


y* 


= Qs 


x,y 


Ey<x0> 


(3) 


The  matrix  Qs  has  elements  which  are  only  functions  of  S^,  while  E^x^) 
is  a linear  ved^or  representing  boundary  conditions  between the  four 
switching  states. 


6 


Solution  in  the  Time  Domain  for  the  State  Variables 


vy(t)  = *y,0  * eXp  ^Sx0  * t)  + ^y,  2 * ®*P  (sxl  * t)  + 

*y,3  * <Sx,2  * t)  Sln  (Sx,3  * * + V (4) 

with 

av  * ^ <SX,3/(SX,2  * A.2»  (5) 

The  subscript  y in  Equations  (4)  and  (5)  defines  the  state  variable  and  the 
coefficients  are  derived  from  Ayx*  The  complete  circuit  analysis  is 
presented  in  Appendix  A and  the  computer  program  for  calculating  the  time 
domain  solutions  is  presented  in  Appendix  B. 

APPLICATION  OF  ANALYSIS 

Steady  State  Conditions 


The  solution  of  the  steady  state  condition  requires  the  cyclic  perio- 
dicity of  the  voltage  waveforms  as  described  by  Equation  (4).  As  shown  in 
Figure  2,  the  sequence  of  the  four  switching  periods  is  fixed.  At  the  end 
of  the  fourth  switching  period,  the  system  must  return  to  the  values  at  the 
beginning  of  the  first  switching  period.  The  mathematical  expression  for 
this  condition  requires  solution  of  a set  of  transcendental  equations.  The 
complexity  of  these  equations  does  not  permit  an  explicit  solution  in  the 
form  of  Equation  (4)  by  defining  the  values  of  Xyx  for  steady  state  as  a 
function  of  the  circuit  parameters.  The  criterion  for  the  existence  of  a 
stable,  steady  state  solution  can,  however,  be  tested  directly  by  execution 
of  the  computer  program  for  any  given  set  of  circuit  input  parameters  C, 

Cq,  L,  R^,  and  M.  A set  of  initial  boundary  conditions  must  be  chosen  for 
the  start.  If  after  a reasonable  length  of  time,  transient  conditions 
cease  and  the  computer  output  shows  a stable  cyclic  wave  train,  the  stable 
steady  state  condition  can  be  directly  read  from  the  specific  features  of 
this  wave  train.  The  criterion  for  a single  valued,  steady  state  condition 
which  is  correlated  with  a set  of  input  circuit  parameters  is  given  by  the 
condition  that  the  steady  state  waveform  is  independent  of  the  initial 
boundary  conditions  selected. 

Mode  1,  Continuous  Inductor  Current  Mode 

An  example  of  the  evaluation  of  a steady  state  condition  for  a given 
set  of  circuit  parameters  and  a chosen  constant  value  Z2  of  the  discrete 
diode  conduction  time  is  shown  in  Figures  3 to  5»  Figure  3 shows  the  first 
six  cycles  after  initial  start-up  condition  with  circuit  input  parameters 
(C,  Cq,  L,  M)  and  operational  parameters  (R^,  Z2)  indicated.  The  wave 
shapes  for  the  normalized  values  of  voltage  Vr  on  the  series  capacitor, 
the  voltage  Vr  on  the  inductor,  and  the  current  J through  the  inductor  are 
plotted  together  with  voltage  VT  on  the  primary  winding  of  the  transformer. 
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Also  plotted  is  the  DC-output  voltage  Vq.  The  scale  factors  for  Vj  and  Vq 
have  been  chosen  such  that  their  ratio  is  identical  with  the  winding 
ratio  M on  the  transformer.  It  is  seen  that  V>p  is  identical  with  the  trace 
of  Vq  in  the  positive  half-cycle  of  J and  Vj  = -M  • Vq  in  the  negative  half 
cycle.  Negative  values  of  VT  are  associated  with  negative  values  of  J. 

The  initial  set  of  boundary  conditions  chosen  for  the  beginning  of  the 
first  conduction  period  of  transistor  TR1  Eire:  = E/2,  J = 0, 


Vrj,  = 0.  This  means  that  the  input  voltEige  is  equally  divided  among  the 
2 series  capacitors  Cl  and  C2  (see  Figure  1).  When  transistor  TR1  is 
turned  on  at  point  Tq,  the  output  filter  capacitor  Co  is  in  a discharged 
condition  so  that  VT(0)  = 0.  The  current  J charges  capacitor  Cl,  while 
capacitor  C2  is  being  dischEirged.  Curve  Vq2  is  not  shown  on  the  plot  but 
is  antisymmetric  to  with  respect  to  a line  drawn  parallel  to  the  time 
axis  with  the  value  E/2.  The  current  J is  almost  sinusoidal,  reaches  a 
maximum,  and  then  swings  toward  zero.  Note  that  the  plotted  interval  of 
points  has  been  reduced  by  a factor  of  10  as  the  current  approaches  zero. 
This  measure  provides  higher  accuracy  for  the  time  limitation  criterion  of 
the  first  transistor  conduction  period  which  ends  when  J reaches  zero.  Ihe 
normalized  inductor  voltage  Vl/E  starts  with  a value  of  0.5  and  swings  to  a 
negative  value  of  approximately  -0.5  at  the  end  of  the  first  transistor 
conduction  period  (point  Tl).  The  voltEiges  Vj  and  Vq  increase  very  slowly, 
due  to  the  fact  that  the  output  filter  capacitor  is  about  10  times  larger 
than  the  resoimnce  capacitor  and  its  equivEilent  value  reflected  into  the 
primary  circuit  is  increased  by  the  factor  l/l&  =25.  This  relationship, 
which  has  been  selected  to  reduce  the  ripple  voltage  in  the  output,  is  very 
important  for  understanding  the  nature  of  start-up  transient  wave  shapes 
during  the  first  few  cycles  which  show  insignificant  dependence  on  the 
value  of  the  load  resistance  R^.  The  conduction  period  of  diode  D1  begins, 
at  point  Tl,  with  a negative  rising  current  J and  ends  when  this  current 
reaches  its  maximum  value  at  T2  Eifter  a time  equivalent  to  Z2  = 0.5.  Kie 
inductor  voltage  is  seen  as  swinging  back  to  zero  Euid  V^/E  approaches  1. 

At  the  beginning  of  the  diode  conduction,  Vj  jumps  to  a negative  value  and 
keeps  rising  in  absolute  value.  At  point  T2,  transistor  TR2  is  turned  on. 
Consulting  the  equivalent  circuit  in  Figure  2,  we  have  a situation  where 
C2  is  dischEirged  (because  Vc^=  E)  and  since  Vm  is  still  very  small , the 
bulk  of  the  input  voltage  E appears  on  the  inductor  as  a negative  voltage 
Vl/E  = -1  as  c£in  be  seen  in  the  graph.  During  the  TR2  conduction  period, 

Vl  approaching  zero  causes  a further  increase  of  the  current  in  the 
negative  direction.  Capacitor  C2  is  rapidly  charged  by  this  current  to  a 
positive  vEilue,  while  capacitor  Cl  dischsirges.  When  the  current  reaches 
its  negative  maximum,  Vcl  and  VL  are  approximately  0 and  Vq2=  E.  As  the 
current  swings  back  to  zero,  Vc2  rises  further  and  reaches  a value  of 
approximately  = 2E  at  the  end  of  the  conduction  period  of  transistor  TR2 


at  point  T3.  Since  Vc2  is  not  shown  on  the  computer  graph,  Figure  3>  its 
value  is  inferred  from  the  condition  Vq2  = E - Vqj_.  With  = -1.1  E, 

Vq2  = 2.1  E.  When  conduction  through  transistor  TR2  ceases  at  point  T3, 
diode  D2  starts  conducting  the  now  positive  current  for  a period  equivalent 
to  Z2.  Note  the  sharp  rise  of  the  diode  current  to  its  maximum  value  at  the 
start  of  this  period  which  is  caused  by  a rapid  discharge  of  capacitor  C2 
from  the  VEilue  of  2. IE  to  approximately  1.1E.  At  point  TU,  transistor  TR1 
is  turned  on  again  and  the  first  cycle  is  completed.  Note  that  the  state 
of  the  circuit  has  now  significantly  changed  as  compared  to  the  initial 
boundary  conditions  at  point  Tq.  Capacitor  Cl  shows  a slight  negative 
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charge  VC1  * -0.1E.  The  initial  current  of  the  next  transistor  conduction 
starts  with  an  appreciable  positive  value  (J  = +0.084  E)  and,  therefore, 
reaches  a much  higher  maximum  value.  The  pr#  cess  of  buildup  of  stored 
energy  in  the  capacitor  and  inductor  continues  until  about  the  end  of 
cycle  3 (Zl  = 12)  where  Vq  reaches  a maximum  value  of  +2.62E.  At  this 
point,  Vq  has  reached  a value  of  O.386E  and  the  energy  dissipation  in  the 
load  resistor  ((O.38  E)2/l00)  is  no  longer  negligible.  The  trend  in  the 
following  cycles  is  now  reversed,  more  stored  transient  energy  is  trans- 
ferred to  the  load  and  dissipated.  Peak  currents  and  voltages  begin  to 
fall.  Figure  4 shows  a larger  portion  of  the  continuing  wave  train  com- 
pressed into  the  time  scale.  It  is  seen  that  after  about  25  cycles  a 
steady  state  condition  is  approached,  which  can  be  best  judged  by  the  slow 
rising  value  of  Vq.  Figure  5 shows  enlarged  values  of  cycles  Zl  = 300~324 
from  which  the  following  results  of  the  steady  state  condition  are  deduced: 

Z2  = 0.5;  ^ = 100  ft  ; zi  = 320 

Tl/2  = 23*T7  /As»  f = 21,03  ^ 

Jmax  /E  = °*°96;  Jav/E  = 0,058;  Jmax/Jav  = 1,644 

VCmax/E  = 2*0TJ  VCpp/E  = 3*14 
VTav/E  = 0.201;  Vq/E  = 1.005 

A second  case  for  a larger  diode  conduction  period  Z2  = 0.75  Is  shown  in 
Figures  6 and  7.  It  is  seen  from  Figure  6 that  the  initial  transient 
maxima  for  Vq  and  J are  now  much  smaller  than  those  for  Z2  = 0.5.  Figure  7 
shows  the  waveforms  under  steady  state  conditions  for  Zl  = 100  - 120  with 
the  following  results: 


Z2  = 0.75;  Rl  = 100  XI  ; Zl  = 120 

Tl/2  = 29,4  Ms;  f = l6*99  kHz 

Jmax/E  = 0,058;  Jav/E  = 0,034>  Jmax/Jav  = 1,663 

W13  = 1‘645;  vcpp/E  = 2,29 

VTav/E  = 0,139;  VE  = 0,695 

A similar  computer  plot  is  shown  for  Z2  = 0.875>  Zl  = 0 to  l6  in 
Figure  8.  Comparing  Figures  3,  6,  and  8 it  is  seen  that  the  start-up 
transients  diminish  when  Z2  increases.  Steady  state  condition  is  reached 
after  about  25  cycles.  Figure  9 shows  the  wave  shapes  for  this  condition 
and  indicates  the  following  parameters : 


Z2  = 0.875;  Rl  = 100  & ; Zl  = 104 

ti/2  = 32,905  Ms;  f = 15,195  ^ 

Jmax/E  = 0,0486;  JJE  = °*02T9;  Jmax/Jav 

V /E  = 1.5443;  V /E  = 2.0886 

Cmax  Cpp 

V /E  = 0.1124;  VjE  = 0.5623 

Tav  O' 


1.740 
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Similar  computer  plots  were  obtained,  for  values  of  Z2  = 0,  0.25,  0.375, 
and  1.  A summary  of  these  results  is  presented  in  Table  1 and  as  a graph 
in  Figure  10.  These  results  are  for  a single  value  of  load  resistance 
(R.  = 100  ft ) which  was  considered  nominal  load.  The  curve  labeled  Vq/E 
in  Figure  10  clearly  shows  the  range  of  voltage  control  Vq/E  to  be  2.5  to 
0.5  for  a range  of  dwell  time  0<Z2<1.  It  is  also  seen  that  the  maximum 
capacitor  voltage  increases  approximately  by  a factor  2 over  this  range, 
while  the  ratio  Jjj/Jav  varies  only  slowly  with  Z2.  In  contrast,  the  maxi- 
mum transient  current  to  steady  state  average  current  ratio  shows  a drastic 
increase  as  Z2  decreases  to  zero  (see  curve  labeled  J^/ JflV ) . This 
behavior  raises  interesting  questions  on  the  dynamic  response  of  this  con- 
verter and  on  the  question  of  how  to  design  a control  system  which  avoids 
excessive  component  stresses.  While  the  basic  tools  for  such  a study  are 
provided  in  this  report,  the  study  of  control  concepts  and  design  optimi- 
zation is  beyond  the  scope  of  this  report. 

Mode  2,  Zero  Diode  Current  Mode 

As  can  be  seen  from  Figure  2 (Case  2 and  Case  4),  whenever  E + Vm>VQ, 
the  diode  (D1  or  D2)  is  negatively  biased  and  current  cannot  flow  in  the 
Inductor  L.  This  condition  occurs  particularly  at  low  load-  r-  , 

initial  start-up  transient  condition.  Figure  11  shows  a computer  graph 
indicating  a step  change  in  the  load  value  from  the  nominal  RL  = 100  ft  to 
= 1000  ft  at  the  end  of  switch  period  Z1  = 210.  All  other  parameters 
are  the  same  as  for  Figure  9 with  the  exception  of  the  time  scale.  Immedi- 
ately after  this  step  change  of  RL,  the  transformer  voltage  Vrp/E  and  the 
output  voltage  Vq/E  start  to  rise.  During  the  diode  conduction  period, 
the  peak  current  value  of  j/E  gradually  decreases,  while  at  the  same  time, 
each  succeeding  peak  value  of  j/E  in  the  transistor  conduction  period  rises. 
Figure  12  shows  a continuation  of  the  graph  of  Figure  11.  At  Z1  = 346,  at 
the  end  of  the  TR2  conduction  period,  1 + Vm/E  exceeds  the  capacitor 
voltage  Vq/E  and  Diode  D2  is  reverse  biased  at  the  beginning  of  Z1  = 3^7. 
This  means  that  the  diode  current  is  zero.  The  transistor  peak  current 
has  reached  its  maximum  at  Z1  = 3^  and  begins  to  fall  in  the  succeeding 
switch  periods.  Since  the  average  value  of  the  transistor  current  is  still 
larger  than  the  equivalent  load  current  through  RL'  the  output  filter 
capacitor  Cq  still  receives  an  appreciable  charging  current  during  the 
transistor  conduction  period  and,  as  a result,  the  values  Vq/ E and  Vt/E 
continue  to  increase.  The  rapid  feedback  between  transistor  current  and 
voltage  on  the  resonance  capacitor  is  seen  between  switching  cycles 
Z1  = 346  to  Z1  = 378.  At  the  next  transistor  switching  period  (Zl  = 382), 
the  transistor  is  reverse  biased  by  the  condition  1 + Vq/E<Vt/e  and 
cannot  conduct.  A new  mode,  Mode  3,  starts. 

Mode  3,  Zero  Transistor,  Zero  Diode  Current  Mode 

This  mode  continues  and  is  characterized  by  an  effective  decoupling  of 
the  primary  input  and  secondary  output  stage.  It  can  only  be  changed  by 
the  slow  decay  of  the  voltage  on  the  filter  capacitor  VQ  through  discharge 
into  the  load  resistor  R^.  Therefore,  it  is  clearly  seen  that  a low  load 
current  (high  value  of  R,)  will  prolong  this  mode.  The  continued  wave 
train,  shown  in  Figure  I3,  indicates  that  at  Zl  = 398  a small  transistor 
current  again  starts  to  flow  when  Vt/E  has  reduced  to  a value  of  0.5. 
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Table  1.  Summary  of  Normalized  Transient  and  Steady-State  Parameters 
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Figure  11.  Transient  Wave  Shapes  for  Z2=0.S75,  Step  Change  of  R- 
from  100-1000  £l  at  Zl=210  L 


Subsequently,  this  current  increases  and  slows  down  the  discharge  of  the 
output  filter  capacitor.  When  the  average  value  of  the  transistor  current 
again  exceeds  the  load  current,  the  output  voltage  rises  again  and  the  slow 
oscillation  of  the  system  between  Mode  2 and  Mode  3 may  continue.  However, 
the  oscillation  will  be  finally  damped  out  as  long  as  is  not  infinite. 

The  curves  in  Figure  13  show  that  for  the  selected  value  of  Z2  = O.875  and 
Rp  = 1000  SI  the  converter  approaches  a stable  state  condition  in  Mode  2. 

CONCLUSIONS 

The  transient  analysis  of  the  series  resonance  type  DC  to  DC  switching 
converter  and  the  implementation  by  a computer  program  have  been  found  to 
be  valuable  tools  for  simulating  the  complex  waveform  train  in  the  various 
switching  modes  of  this  converter.  The  free  running  feature  of  this  con- 
verter is  simulated  by  the  restrictions  of  unipolar  current  flow  in  the 
switching  transistors  and  antiparallel  diodes.  The  only  control  parameter 
is  the  diode  conduction  time  (Z2)  which  separates  the  "on  state"  of  the  two 
switching  transistors.  This  leads  to  3 operational  modes  of  the  converter, 
distinguished  by 

Mode  1 inductor  current  conduction  during  both  transistor 
and  diode  conduction  periods; 

Mode  2 inductor  current  conduction  only  during  transistor 
conduction  period; 

Mode  3 zero  inductor  current  during  both  transistor  and 
diode  conduction  periods. 

Mode  3 is  an  intermittent  state  only  but  would  apply  for  the  border  case 
of  an  idealized  converter  (no  component  losses)  under  a no  load  condition, 
which  would  be  an  undefined  case  with  no  steady  state  condition.  For 
finite  values  of  the  load  resistance,  the  converter  has  stable  steady 
state  conditions  in  either  Mode  1 or  Mode  2 depending  on  whether  a heavy 
or  a light  Load  is  applied.  The  transient  Mode  3 appears  Under  a light 
load  condition  with  an  increasing  duration  as  the  no  load  condition  is 
approached. 

For  each  value  of  diode  conduction  period  0<Z2<1,  a single  set  of 
normalized  operational  parameters  can  be  shown  to  exist  in  the  steady  state 
mode.  These  parameters  define  the  available  range  of  the  ratio  of  input 
voltage/ output  voltage,  electrical  stresses  in  the  power  circuit  components, 
and  operating  switching  frequency.  Since  the  average  value  of  the  primary 
voltage  on  the  transformer  cannot  exceed  a value  of  0.5E  (half  the  DC- 
input  voltage),  in  the  steady  state  condition,  the  converter  may  be 
regarded  as  a buck  type  switching  converter  with  a limited  maximum  output/ 
input  voltage  ratio  of  0.5. 

The  computer  program  also  yields  the  maximum  transient  voltage  stresses 
as  a function  of  Z2  and  allows  determination  of  the  settling  time  toward 
steady  state  conditions. 
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Transient  stresses  are  shown  to  be  a maximum  for  Z2  = 0 and  will 
almost  completely  subside  for  Z2  = 1.  Therefore,  control  of  the  parameter 
Z2  can  be  used  effectively  to  suppress  transient  stresses  for  any  opera- 
tional condition  except  excessive  input  voltage  that  would  have  to  be 
contained  by  an  appropriate  input  transient  suppression  circuit.  The 
analysis  and  computer  simulation  considers  lossless  power  train  components, 
therefore,  it  gives  no  information  regarding  the  operational  changes  of 
converter  efficiency. 

Two  considerations  for  component  losses  may  be  made  to  assess  devia- 
tions from  the  ideal  case: 

1.  Current  independent  losses,  such  as  transformer  core  losses 
and  leakage  currents  in  the  output  rectifier  and  filter  capacitor,  may  be 
regarded  as  an  additional  shunt  loss  to  the  output  load.  Therefore,  its 
prime  operational  influence  will  be  to  remove  the  zero  output  load 
instability  and  present  the  no  load  case  as  a Mode  2 operational  condition. 

2.  Current  dependent  losses  like  copper  losses  in  the  inductor 
and  transformer  can  be  simulated  as  a lumped  resistance  in  series  with 
the  output  load.  In  general,  if  these  losses  do  not  exceed  10$  of  the 
power  rating  of  the  converter,  their  influence  on  the  normalized  parameters 
is  minor  and  does  not  cause  significant  deviation  of  the  operational  mode 
of  the  converter  as  presented  in  the  analysis  of  the  ideal,  lossless  con- 
verter. Therefore,  in  most  design  cases,  the  ideal  component-converter 
analysis  will  provide  a good  representation  for  converter  operation. 

Power  factor  losses  will  increase  with  the  ratio  Jm/Jav  as  Z2  increases, 
but  as  this  ratio  changes  only  by  20$  for  0<Z2<1  there  is  no  dramatic 
effect  on  efficiency  expected  over  a wide  voltage  control  range. 

The  computer  program  is  very  versatile  and  permits  parametric  studies 
for  design  optimization  by  varying  either  the  circuit  parameters  or  the 
switching  control  parameter. 

RECOMMENDATIONS 

The  application  of  the  computer  program  for  analysis  of  control  loop 
design  and  optimization  for  the  DC/DC  switching  converter  is  recommended. 

It  is  further  recommended  that  a study  be  conducted  to  determine  the 
advisability  of  replacing  the  ideal  lumped  parameter  presentation  of  power 
components  and  switches  by  more  refined  component  models.  A trade-off 
determination  should  be  made  between  the  increased  complexity  of  the 
computer  program  and  the  necessity  for  more  accurate  and  realistic  results. 
The  latter  study  will  be  a necessity  when  further  increases  in  switching 
frequencies  beyond  the  100  kHz  range  are  anticipated. 
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Circuit  Analysis 


The  equivalent  circuit,  which  is  derived  from  the  simplified  power 
switching  circuit  of  the  series  capacitor  resonance  DC  to  DC  converter,  is 
shown  in  Figure  A-l.  It  represents  a DC  input  voltage  E,  a switch  S,  which 
can  be  opened  or  closed  in  accordance  with  the  operational  switching  require- 
ments of  the  converter,  a voltage  generator  v^,  which  replaces  the  passive 
circuit  configuration  of  the  output  circuit  as  reflected  into  the  primary  of 
the  converter  transformer,  the  resonance  inductor  L,  and  the  resonance  capac- 
itor C.  Figure  A-2  shows  the  equivalent  circuit  for  vrp  on  which  the 
restriction:  1 


vT  = sign  (i)  . v0  • M 

(A-l) 

is  imposed. 

i • . . . . transformer  primary  current 
v0  ....  DC  output  voltage 

M ....  W1/W2  turns  ratio  of  converter  transformer 

Hie  following  set  of  simultaneous  equations  describe  the 

state  of  the 

circuit  at  a particular  time  instant  t. 
E - vT  - vL  - vc  = 0 

(A-2) 

vL  = L • di/dt 

(A-3) 

vc  = (l/c)  * J*i  • dt 

(A-l*) 

v?  = sign  (i)  • vQ  • M 

(A-5) 

v0  = vo(°)  + (VCq)  • • dt 

(A-6) 

icQ  = abs  " vo/rl 

(A-7) 

Differentiation  of  (A-6)  together  with  (A-7)  yields 
dvQ/dt  = (1/CQ)  • [abs  (i/M)  - Vq/r] 

( A— 8 ) 

and 

i = M • C0  [dvQ/dt  + vq/(Rl  • C0)j  sign  (i) 

(A-9) 
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Substitution  of  (A-9)  into  (A-3)  yields: 

vL  = L • M • C0  [d2Vo/dt2  + dv0/dt  (1/RL  • C0)] 
Substitution  of  (A-5)  into  (A-10)  yields: 

= L • CQ  d^/dt2  + (L/R^)  • dvT/dt 

Differentiating  (A-4)  yields: 
dvc/dt  = i/C 

Integration  of  (A-3)  yields: 

L • i = • dt  + vl(0) 


where 


1 -£  /''Ldt*  ‘l*0* 
i (0)  = vl(0)/L 


Combination  of  (A-12)  + (A-15)  yields 


dvc/dt  = i • [l  r /vL dt  + iL(0)] 


and  differentiation  of  (A-l6)  yields 

d2vc/dt2  = vL  • 1/(L  • C)  (A- 17) 

Equations  (A-2),  (A-ll),  and  (A-17)  provide  a system  of  linear  differential 
equations  with  the  circuit  voltages  vT,  vL,  and  vc  as  variables.  These  are 
chosen  as  the  components  of  the  state  variable. 

V(V  V V °r  v(vi’  V V 


V1  = V V2  = V V3  = VC 

The  set  of  second  order  linear  differential  equations  (A-2),  (A-ll),  and 
(A-17)  transferred  into  3-space  (Laplace  transformation)  yields  the  vector 
equation: 


where : 


V(y)  • S(x,y)  = F(x) 


S(x,y)  state  function  matrix  for  S-space 

F(x)  vector  defined  by  boundary  conditions  at  time-zero 


! 
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It  follows  from  equations  (A-2),  (A-ll),  and  (A-17)  that: 


S(x,y)  =(  S + S/Tq  -l/(T  • Tx)  0 

\ 2 2, 

\ 0 -W  S 

with  Tq  = • CQ;  - L/R^ja)2  = l/(L  • C)  and: 

(e/s  \ 

Vx  (0)(S  + 1/T0)  + vi  (0)  J 
s • v,(o)  + v'  (0)/ 


(A-19) 


(A-20) 


where:  V-|_(0)  and  V_(0)  initial  values  of  the  state  function  component  vectors 

and  V^(0)  and  V (o)  initial  values  of  the  derivative  of  the  state  function 
^ component  vectors  and  V^. 

The  Eigenvalues  or  roots  of  the  characteristic  equation  follow  from  setting 


Determinant  S(x,y)  = 0 which  yields: 

S4  + S3/Tq  + E?(1/Tq  • Tx)  + ft)2)  - S ‘ft)2/T0  = 0 (A-21) 

Equation  (A-21)  is  of  fourth  order  with  one  root  SQ  = 0.  The  remaining 
cubic  equation 

S3  + S2/Tq  + S(1/(T0*  T^)  +ft)2)  + ft)2/T0  = 0 (A-22) 


is  characteristic  for  a dumped  oscillation  with  the  angular  frequency  a) 
and  damping  constants  depending  on  Tq  and  T-^. 


For  the  oscillatory  mode,  the  other  roots  (Eigenvalues)  of 
equation  (A-22),  labeled  S^,  So,  S^  can  be  written  as:  S,;  Sp  = 
= su  - js5. 


Si*  + J 


From  equation  (A-l8)  it  follows 

V(y,S)  = F(x)  • S_1(x,y)  (A-23) 

which  can  be  written  in  the  form 


v(y,s) 


x.L°iXl  ♦ ZLhxl  * 

S - SO  S - SL 


X(2,y)»s  4 X(3,y) 

(s  - su)2  + s52 


(A-24) 
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Equation  (A-24)  yields  the  general  form  of  the  components  of  the  state 
variable  in  the  time  domain  as 


V(y,t)  = X(0,y)  + X(l,y)  • exp  (S  • t)  + X(2,y)  i 

S5 

[(S.  + X(3,y)/X(2,y))2  + S2]  1/2  • exp  (S.  • t)  • 


sin  • ( • t + a) 


where  a = arct  • (S5/(S^  + X(3,y)/X(2,y) ) 


(A-25) 


(A-26) 


Derivation  of  Constants  X(n,y) 

Ihe  constants  X(n,y)  are  defined  by  equations  (A-23)  and  (A-24). 
Equation  (A-2U)  can  be  written  in  the  general  form: 


, , A^(s) 

V(y,s). 


(A-27) 


where  A3(s)  is  a third  order  polynom  of  s and  b\s)  is  a fourth  order 
polynom  of  s with  the  value  B^(s)  = Det.  S(x,y)  and  SQ,  S1?  S^,  and  Sc  are 
roots  or  eigenvalues  of  S(x,y)  (Equation  (A-21)). 

On  the  other  hand,  (A-23)  can  be  written  as: 


V(y,s)  = Det.  Sp^^x^/Det.  S(x,y) 


(A-23) 


where  is  formed  by  replacing  the  y-column  of  S(x,y)  by  the  column 

F(x).  It  follows  from  (A-27)  and  (A-28)  that 


A3(s)  = Det.  Sp(x)(x»y) 

Since  (A-29)  must  have  equal  coefficients  of  sn  on  both  sides  of  the 
equation,  it  follows  from  (A-24)  and  (A-29)  that: 


(A-29) 


X(0,y)  (s-sL)  [ (s-s4)2  + s52]  + X(l,y)  • (s-sQ)  [(s-s^)2  + s52]  + 

X(2,y)  • (s-s0)  (s-s1)  s + X(3,y)  (s-s0)  (s-s-^  = Det.  Sp(x)(x>y)  (A-30) 
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Applying  the  aforsnentioned  condition  of  equal  coefficients  of  sn  to 
equation  (A-30)  yield  a set  of  k linear  equations  in  X(n,y)  as  follows: 


r 


X(0,y) 


+ X(l,y) 

(-s!  " 2s4) 

+ X(l,y)  • 

(-s0  - 2sk) 

[(s42  + s52)  + 2sis4)] 

+ X(l,y)  • 

[(S42  + s52)  + 2s0sJ 

[c%2  * s52>  • (-v  ] 

+ X(l,y)  • 

[(s42  + s52)  * ("so)] 

+ X(2,y) 
+ X(2,y) 

+ X(2,y) 
+ 0 


+ 0 

(-sl  - sQ)  + X(3,y) 


- h 
= E2 


(sQ  • s ±)  + X(3,M  (-30  - sx)  = ^ 
+ X(3,4)  • s0s1  = 


or  in  matrix  form: 

X(n,y)  • Sn(x,y)  = E(n,y) 
where  the  matrix  S^Xjy)  follows  from  (A- 31)  as: 


(A-31) 


(A-32) 


■ 1 
! 


r 1 

1 

1 

0 

1 

“(s1  + z\) 

“(8o  + 2s0 

-(31  * 8o> 

1 

1 

Sn(x,y)  = 

s42  + s52  + 2sis4 

3k  + s52  + 2s0s4 

S0S1 

-(s0+Sl) 

(A-33) 

l] 

fl 

(^k2  + s52)  • (-3l) 

(31|2  + s52)  (-8o) 

0 

8081 

J 

• 

and  the  matrix  E(n,y)  is  given  by  the  coefficients  of  sn  for 
n = 3,  2,  1,  0 when  Det.  is  evaluated  as  a polynomial  of  sn 
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vx(0) 

Vl(°)/T0  + v3  (0) 

E/(T0  • Tl)+  2v1(0)  - v3(0)/(TqT1) 
v1(0)  • U)2/T0  + <11^3  (°)  “ v3  (0)/(T0  • Tx) 


E - vi(0)  - v3(0) 

E/T0  - (vx(0)  + v3(0))/Tq  - (0)  - v^  (0) 


V (0)/Tq 


v3(0) 

v3(0)/Tq  + v’  (0) 


W2  E - v1(o) j + v3(0)/(Tq  • Tx)  + v3  (0)/T0 
u>2’  E/Tq  - v^OJ/Tq  - (0)]  + (Oj/T^) 


from  equations  (A-32),  (A-33)>  and  (A-3M>  it  follows 

X(n,y)  = E(n,y)  • S^1  (x,y)  (A-3" 

Using  equations  (A-25),  (A-26),  and  (A-35),  the  state  function  V(y,t)  is 
defined  in  the  time  domain  in  terms  of  circuit  parameters  which  enter  into 
the  expressions  for  TQ  = RLCQj  T,  = Rj/L  andw=  l/(LC)1/2  and  in  terms  of 
initial  boundary  conditions  for  the  transformer  primary  voltage  v1(0)  and 
the  capacitor  voltage  v^O)  and  their  derivatives  v|  (0)  and  v^  (0)  which 
are  contained  in  matrixJE(n,y) . 
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APPENDIX  B 


Computer  Program  to  Analyze  and  Plot 
Waveforms  in  the  Time  Domain 


INTRODUCTION 

The  computer  program  is  written  in  Basic  4 language  as  implemented  on 
a time-shared  system  (TSS/3)  using  a PDP-8/l  Minicomputer.  The  organization 
of  the  program  closely  follows  the  mathematical  modeling  analysis  of  the 
switching  type  converter  as  described  in  Appendix  A. 

The  program  is  subdivided  in  4 segments,  which  are  stored  on  the  system 
device  (RF-08  Disk)  in  the  user  protected  area  and  called  into  core  for 
operation  as  required  by  the  program.  An  independent  data  record  file  has 
been  established  for  the  program  -which  handles  input  and  output  of  data  in 
the  4 chained  segnents  of  the  program.  Small  supplementary  utility  programs 
were  written  to  print  the  content  of  the  data  file  and  make  changes  of 
certain  circuit  parameters  for  the  purpose  of  studying  transient  response 
phenomena. 

The  program  senses  for  negative  bias  condition  of  the  transistors  and 
diodes.  When  this  condition  occurs,  the  inductor  current  is  set  to  zero 
and  the  converter  operates  in  a first-order  mode  defined  by  the  discharge 
of  the  output  filter  capacitor  into  the  load  resistance. 

PROGRAM  ORGANIZATION 

The  first  program  segment  has  the  file  name  EKPPO.  It  provides  the 
initial  input  to  the  data  file  which  contains  the  circuit  parameters  C,  CO, 
L,  Rl,  and  M and  the  specified  discrete  value  Z2  for  the  diode  conduction 
time.  Additional  parameters  like  Zl,  the  counting  variable  for  the  string 
of  consecutive  switching  periods,  the  constant  Z5  which  sets  an  upper  limit 
for  the  value  of  Zl  in  the  execution  of  the  program,  and  scaling  parameters 
for  the  graphic  plotter  are  also  defined  in  this  first  program  segment. 

The  major  function  of  segment  EKPPO  is  to  define  and  solve  the  charac- 
teristic equation  of  the  circuit  and  provide  the  roots  of  the  characteristic 
equation  as  input  to  the  data  file  for  further  processing  in  the  follow-on 
program  segments  (see  line  statements  l8l  through  193)* 

Ihe  second  program  segment  has  the  file  name  EKPP1.  Its  main  function 
is  to  define  the  boundary  conditions  for  each  of  the  4 switching  states  in 
terms  of  either  arbitrary  initial  circuit  conditions,  or  as  determined  from 
the  state  of  the  circuit  at  the  end  of  the  previous  switching  state.  The 
matrix  elements  of  the  state  function  equation  are  calculated  and  stored 
in  the  data  file. 

The  third  program  segnent  has  the  file  name  LINEQ1.  This  is  a sub- 
program to  invert  the  matrix  of  the  state  function  and  calculate  the  set 
of  constants  to  be  used  with  the  time  domain  solutions  of  the  state 
variables. 
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The  fourth  and  final  program  segment  has  the  file  name  EKPP3.  In  this 
part  of  the  program,  the  time  dependent  values  of  the  state  variables  are 
calculated  and  printed,  and/or  plotted  on  a graphic  terminal.  Hie  values 
of  the  state  variables  at  the  end  of  the  switching  period  are  stored  in  the 
data  file  for  processing  of  boundary  conditions  for  the  subsequent  switching 
period.  In  case  of  continuous  plotting,  this  program  segment  is  chained  to 
program  segment  EKPPl.  The  program  loops  through  the  segments  EKPP1,  LINEQ1, 
EKFP3,  EKPPl,  etc.,  in  a continuous  fashion  until  the  specified  number  of 
switching  periods  Z5  is  reached.  The  program  also  contains  statements  to 
limit  the  transistor  and  diode  switching  periods  whenever  the  primary  loop 
current  reaches  zero,  thereby  modeling  the  conditions  for  unipolar  conduc- 
tion of  transistors  and  diodes.  Thus,  the  conduction  time  is  independent 
of  the  specified  range  of  the  time  variable  which  does  not  exceed  the  hwl f 
period  time  of  the  resonant  frequency  but  it  can  be  much  shorter  in  the  case 
of  specific  operational  modes  of  the  converter. 


Two  utility  programs  have  been  written  which  provide  for  checking  and 
new  input  to  the  data  file,  independent  of  the  program  execution. 

Program  segment  "EKD1"  prints  the  content  of  the  data  file  as  output 
on  the  teletype.  Program  segment  "EKD2"  is  used  to  change  the  input 
parameters  of  the  circuit  (C,  CO,  L,  Rl,  and  M)  and  the  initial  boundary 
conditions  and  plotting  range,  as  defined  by  data  lines  numbers  6l  through 

68. 


If  the  program  is  halted  in  any  switching  period,  the  combined  use  of 
EKD1  and  EKD2  can  facilitate  checking  of  boundary  relations  and  input  of 
step  functions  for  the  purpose  of  studying  transient  behavior. 

COMMENTS 


A final  note  on  the  usefulness  of  this  computer  modeling  program  is  in 
order.  Safeguards  have  been  provided  to  abort  the  program  if  circuit  param- 
eters are  selected  which  do  not  meet  the  criteria  for  resonance  oscillation 
in  the  primary  circuit. 

In  program  segment  "EKPPO"  there  is  a printout  of  the  critical  load 
resistance  R2.  A value  Rl<  R2  would  provide  overdamping  without  oscilla- 
tions and  an  abortion  of  the  program.  On  the  other  hand,  a value  Rl  » R2 
will  cause  long  transient  periods  when  operational  parameters  are  changed 
and  in  some  specific  cases  will  result  in  modes  during  which  either  the 
diode  or  the  transistor  current  or  both  are  zero  during  a number  of  con- 
secutive switching  periods. 


A study  of  these  conditions  is  of  particular  interest  with  respect  to 
stability  problems  of  the  circuit  but  is  beyond  the  scope  of  this  report. 

Although  the  program  has  been  written  for  a constant  value  of  Z2,  the 
program  can  be  interrupted  at  any  time  to  change  the  value  of  Z2  and  to 
continue  the  study  of  transient  behavior  with  this  new  value  of  72. 
Ultimately,  Z2  could  be  made  subject  to  restriction  of  output  parameters 
(output  voltage,  output  current ) or  limitation  of  stresses  in  electrode 
components  (VQpp)  to  implement  desired  control  functions  with  automatic 
feedback. 
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LISTING  OF  PROGRAM  SEGMENTS 


The  following  is  a list  of  symbols  used  in  the  various  program 


segments: 

A(0,0)  to  A(3,3) 

Matrix  elements  of  state  function  matrix 

B(0)  to  B(23) 

Matrix  elements  of  sub-matrices  of  state  function  matrix 

C 

Equivalent  resonance  capacitance  (C  = Cl  + C2) 

CO 

Output  filter  capacitance 

D (in  EKPPO) 

Criterion  for  solution  of  characteristic  equation 

DO 

Criterion  for  zero  current  approach 

d4(i) 

Array  variable  (I  = 0 to  20) 

D^(o) 

Valuable  to  store  capacitor  voltage  VQ(l) 

D4(l) 

Variable  to  store  inductor  voltage  VQ( 0) 

Dk(2) 

Variable  to  store  inductor  current  JQ(l) 

D4(3) 

Dummy  variable  (not  used) 

D4(lf) 

Variable  to  store  elapsed  time  (T2) 

Di+(5) 

Variable  to  store  capacitor  voltage  V^(l) 

D4(6) 

Variable  to  store  inductor  current  Jq(1) 

D4(7) 

Variable  to  store  inductor  voltage  VQ(0) 

D4(8)  thru  d4(20) 

Dummy  variables  (not  used) 

D5(I) 

Array  variable  (I  = 0 to  7) 

D5(0) 

Variable  to  store  state  function  subscript  Y 

D5(l) 

Variable  to  store  elapsed  switching  time  N1 

D5(2) 

Variable  to  store  capacitor  voltage  VQ(l) 

D5(3) 

Variable  to  store  inductor  current  J^(l) 

D5(U) 

Variable  to  store  transformer  voltage  V^(2) 

D5(5) 

Variable  to  store  switch  count  variable  Z1 

D5(6) 

Variable  to  store  capacitor  voltage  VQ(l) 

D5(7) 

Variable  to  store  inductor  current  JQ(1) 

Array  variable  (I  = 0 to  3) 

Initial  value  of  capacitor  voltage  Vq(1) 

First  derivative  of  capacitor  voltage  at  time  zero 
Initial  value  of  transformer  voltage  VQ(2) 

First  derivative  of  transformer  voltage  at  time  zero 
Array  variable  (I  = 0 to  3) 


Variables  defining  switch  boundary  conditions 


J0(1) 


Array  variable  (I  = 0 to  4) 

Storage  variable  for  capacitance  value  C 
Storage  variable  for  capacitance  value  CQ 
Storage  variable  for  inductance  value  L 
Storage  variable  for  load  resistance 
Storage  variable  for  transformer  winding  ratio  M 
Array  variable  (I  - 0 to  6) 

Storage  variable  for  root  S^_  of  characteristic  equation 

Storage  variable  for  root  S^  of  characteristic  equation 

Storage  variable  for  root  S^  of  characteristic  equation 

Storage  variable  for  operating  frequency 
Storage  variable  for  time  constant  TQ  = • CQ 

Storage  variable  for  time  constant  = R^/L 

Storage  variable  for  angular  frequency  W4  = 1/ Vl»C 
Integer  variable  used  in  array  statements 
Integer  variable  used  in  array  statements 
Variable  describing  inductor  current  normalized  value 


K(I) 

K(0) 

K(l) 

K(2) 

K(3) 

L 

L(I) 

M 

N1 

N2 

P 

Q 

Q2 

03 

Q4 

R1 

R2 

51 

52 

53 
TO 
T1 
T2 
T5 
T6 


Array  variable  (I  * 0 to  3) 

Stores  progressing  elapsed  time  of  each  switching 
period  (Nl) 

Specifies  number  of  switching  periods  Z5 

Full  scale  value  of  time  for  graphic  plot 

Stores  normalized  diode  conduction  time  Z2 

Value  of  inductance  in  series-resonant  circuit 

Dummary  array  variable  (not  used)  (I  = 0 to  b) 

Value  of  transformer  winding  ratio  Wl/W2 

Value  of  switching  period  time 

Time  increment 

Parameter  in  cubic  equation 

Parameter  in  cubic  equation 

Parameter  in  cubic  equation 

Parameter  in  cubic  equation 

Parameter  in  cubic  equation 

Parameter  in  cubic  equation 

Value  of  load  resistance  R^ 

Value  of  critical  load  resistance 

Root  of  characteristic  equation 

Root  of  characteristic  equation  (real  part) 

Root  of  characteristic  equation  (imaginary  part) 

Time  constant  (TO  = R^  • CO) 

Time  constant  (T1  = RjVl) 

Time  variable  (within  switching  period) 

Total  elapsed  time  (at  end  of  switching  period) 
Full  scale  value  of  time  for  graphic  plot 


— ■-  


I 

i 

i 


i ■ 
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T7 


I 

.1 


1 


T(3,3) 


VO(I) 

vo(o) 

V0(1) 

V0(2) 

V0(3) 


W3 

w4 


X$(I) 

X(I) 

x(o) 

X(l) 

X(2) 

x(3) 

Y 

Y1 

Y2 

Z 

Z1 

Z5 


Upper  limit  value  of  time  for  graphic  plot 
Array  variable  (matrix  T(3>  3) = Matrix  A“^(3>3) 
Array  variable  (I  = 0 to  3) 

Variable  inductor  voltage 
Variable  capacitor  voltage 
Variable  transformer  voltage 
Variable  output  voltage 
Effective  angular  frequency 
Resonance  angular  frequency 
String  variable  (I  = 0 to  3) 

Array  variable  (I  = 0 to  3) 

Constant  in  time  dependent  solutions  of  VO 
Constant  in  time  dependent  solutions  of  VO 
Constant  in  time  dependent  solutions  of  VO 
Constant  in  time  dependent  solutions  of  VO 
Index  of  state  function  variable 
Svitch  mode  index 
Switch  mode  index 

Argument  for  Sin  and  Cos  functions 

Counting  variable  for  switch  periods 

Number  of  switch  periods(Zl)  specified  for  plotting 
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PROGRAM  SEGMENT  "EKPPO' 


'•] 

I 


• h t-ASICA 

NEk  OK  01D--0IP 

OLI:  EKOCKAK  NAKE--EKEKH 


: 


n 

EAPY 


LIST 

4 KECOKU  C< 4) ,H<  6) ,Ob< 7)*H(4) 

6 OPEN  8#"PA1A1" 

8 l ET  J =0* 

31  b'K  I NT 

tv  input  r,c0#i  ,ki 

a)  1 FI  G(0)=C\LET  CC  I )*C&M  El  G(2)=IMEI  CC3)=nl\IM  G(4)=p 
4b  LE  T TP=P  »2*K  1 * (C0/P  »2  ) 

*6  LFT  7 I =1  7'P  *2**  1 ) 
bP  Phi  NT  "IP*"  1 T0  J"T  I ="  IT  I 

bb  LET  t 3=1  /(  SPhCl  *C*(C0/k  f2  ) 7<C078  *2*C  ) ) ) 

60  LET  V *=l /SI  K ( L *0) 

61  PnlNT  "V3  = "V.3  J"U="  JM 
6b  Ph  1 M 

66  FhINl  "KEDl'CED  CUPIC  ElUATION  S*3*P*S*t*0  " 

70  l E 7 Psn.1»?)-I/(3*TP»?) 

7b  IET  C =(2*9*(T0*2)*(3*( E 4 f 2 ) - ( h 3 *2  ) ) )/(3*T0) *3 

8b  IET  0=(P/3) *3*(L/?)t? 

87  LEI  ft?  = l /(U4*CP*SC  KC  3)  ) 

88  PKINT  "K2="Jh2 

90  IF  SON(D>«=0  THEN  100 
9b  GOTO  1 I 0 

100  PhJNT  "KESl'l  7 OVEhDAPPED  OSCJLLA1JON  CNANCF  CJnClIl  0 An AP  E I F K b" 
10?  INPUT  C#C0#l ,KI  #P 
10b  GOTO  4b 

110  PH  INI  "OSCILIATOKY  CASE" 

111  PH  INI 

12b  LET  t 3=  - C 72  ♦ SC' n ( D ) 

1 ? 6 I El  C'As-C  /2-SC  h <D> 

127  ] F F 3 <0  THEN  l?9 

128  GOTO  133 

129  I ET  Gl=(-G3)*( I 73) 

I 30  LET  < I =-(  I 

132  COlO  13a 

1 33  LET  G I *G  3 » (173) 

134  IF  64«0  THFN  136 
1 3b  GOTO  I 40 

136  LET  C 2 * ( -( 4)  * ( 1 73 ) 

137  LET  G2«-C2 
139  GOTO  16b 

I 40  LET  C'2*E4»<|/3> 

16b  LET  SI  =(C  I ♦(  2 )-  1 /<  3*  10. ) 

170  LEI  S2*-SI 72-1 /(?• 10  ) 

171  LET  S3  = SCH(3)*<t 1 -(  2)72 

180  PPINl 

181  PKINT  "KESlll  T" 

IP?  PKINT  "THE  HOOTS  OF  THE  CHAHAC 1 Ek I ST  I C t ULA T I ON  AhE  : " 

183  PRINT  "REAL  h 00 7 SI*  "l SI 

184  PKINT  "C08P1F*  rcOOT  S?  ♦ J*  S3  *"  I S2  I" ♦ J*"  I S3 
IPb  PRINT  "OOPPIF.X  HOOT  S2 - J* S3 «" I S? I" - J*" I S3 
I 88  LET  F0*S3 7(2*3. I 4| b93> 

190  PhINT  "OPEKATING  FPE(UFNCY  F 0 •"  I S3/<  2*3 . 1 41  b93  ) l"HF  ft  1 7" 

19?  LET  8 ( 0 ) *S  I M ET  N(|)«S?\LFT  HC2)«S3MFT  H(3)*F0 

193  LET  H(4)*T0MFT  H(b)*tlMFl  H(6)*W* 

194  LET  r>b(P>*PMF1  0b(  b)  *0\H  (0  )*0%8(  I )*0N8(2  )*0N8(  3)  *0 

19b  PKINT  "NUPPFK  OF  SWIICHINI  PFMOPS  TO  PI  01  7b*"l\|NPUl  7b 
196  k ( I ) *7b 

19/  T 7 = 7b*3 . I *1  b93/HC? >\Ph | N|  " 1 7 *" I 1 7 

198  PKINT  "FUllSfAIF  V .1  UF  OF  1IPF  A*lb  I6«"IWNPUT  16*8(21*16 

199  PKINT  "PIOOF  CONU.  -CT  | ON  PF n | OP  7?*"lN|NPl»T  72\k(3)*72 

200  l FT  1*0 
20?  PUT  8,*,J 

210  CHAIN  "FkPP I " • 

2b0  F NO 


36 


j 


wr 


PROGRAM  SEGMENT  "EKPP1" 


1 


1 

I 


01  P 

01  D FhOCnAF  NAKE--EKFFI 
K 

EAPY 

I JM 

b2  HFCOhD  CC4).HC6>.PSC 7>.KC4> 

62  OFEN  P*"DATAI" 

72  LET  F = 2 
P2  CM  P*  b2#  I 

F7  I FT  M*HC2)\LFT  S2*HCI>MF1  S3*HC2>\LET  h2*HC3>\LET  12*HC4> 

PP  IM  TI*WCb>\lFT  k4*H(6)MFl  72*MC3> 

132  KFCOhD  «(3»3)#F(3),l(4) 

136  DIF  MOFiEOF 

137  K OF  I *2  TO  3 
I 3P  ft E A D X SC  I > 

I 42  NFX1  I 

141  Y*P5C2>  \ IF  Y»2  THEN  143 

14?  FhINT  "INITIAL  ROIINDAHY  CONTITIONb  AT  TIKE  ‘ F " 

1 43  7l=nS<5>\|F  7 1*2  THIN  I 4b 

|44  COTO  >46 

|4b  IF  Y=2  THEN  162 

146  73*7l-4*INTC7|/4>\YI*C-l>»lNTC73>\Y2*<-l)»INTC73/2> 

147  IF  Y I = 1 THEN  149 
14P  IF  Y 1 * “ I THFN  1 b2 
»49  IF  Y? - 1 THFN  1 b4 

I b2  IF  Y?*-l  THFN  I b4 

1S1  Ff  I NT  " INFt'T  Fh  h Ok  " MOTO  144 
lb?  IF  Y2*-1  1HFN  I b9 
I b3  IF  Y?*l  THFN  I b9 

Ib4  IF  72*2  THFN  262 

Ibb  Eff’)*1“0bf6)\FC  I > *-P5t 7 >/0C 2 > NEC? >*-DbC 4) 

I b6  EC3 ) = -0b( 7)*C<4 ) t?/C( 1 )*SCN(Db< 7) >*AF SCE(2 ) )/H< 4)\|F  Y*F  1 HEN  172 
1 b7  Fh  I NT  "H  01  IN  PA  HY  CONDITIONS  FOh  ThANb.  CONDUCTION  FEhlOD  7 I ■"  I 7 I 
IbP  FhINT  "FCF)*”IF  < C >#'*F  ( 1 ) *"  I E C I >»,,EC2>*"1EC2)#,*EC  3>***JEC3>NG0T0  I 72 
I b9  FC?  >*-PbC  4>\F  C3  >s-FC?  )/HC  4>NFC2)sDbC2  >NF C I >*2\IF  Y *2  THEN  I 72 


1 6k 

F.t  I NT 

"FOUNDAhY  CONDITIONb  FOh 

DIODE  CONDUCTION  FEhlOD  71* 

"1?  1 

161 

FhINT 

"M  2 > *"  1 F <2)*"E<I  )*"IF<  1 ) 

,"E 

c?  >*•• 

IF  C?  >i 

."E  « 3 ) * 

••i  C 3 ) \C  01  0 

1 6'> 

IF  Y » 

2 THFN  Ibb 

1 63 

FhINT 

"VOLTM?  3\  CAFACITOh  Cl 

■ •1 

1 IKE 

71  hO 

F C2  )«” 

• ' 1 NE  L T 

E 1 » ' 

1 64 

FKINT 

"1ST  DEh I VAT | VE  OF  EC  2) 

AT 

T I v • 

> 0 

E C 1 )*" 

IMNFLI 

EC  1 ) 

l 6b 

Fh  IM 

"VOL  T6(  ‘ f\  ThANbF . CFhlF  ) 

AT 

TIFF 

7 EhO 

E C 2 > ■" 

IMNFL1 

Mr  > 

I 66 

FhINT 

"IbT  DERIVATIVE  OF  E C 2 ) 

A 1 

Iff 

7EhO 

EC  3)«" 

IMNKUI 

E C 3 > 

167  Db< 6>*1 -F <1  ’ b* 7)*-F < 1 >\Db< 4>*-F f 2 >\LOTO  172 

17?  IFT  «r,f)*1MET  AC2,1>«I\LM  AC2#2>«|NLET  AC?,3>*l 

17?  LET  AC  I *2  >*-SI  -2*S?M  ET  A<  | , | ) «-2*b?\L  ET  A C I , 2 > * - b I NLE I A(|,3)«l 

174  IM  A C ? » ? >*S?*?»S3*?*2*bl  • S2  Nt  E I AC  2, I > *S2 »2*b3 *2 

176  IFT  AC?#?>*2MET  AC?,3>«-bl 

|7h  LET  Af  3*2  )*-bl  •<b?*24S3*2)NLET  AC  3,  I ) *F«  NLfc  1 AC3#2>*2 
179  IFT  AC3*3  >*?  NFh  | NT 
I P 4 IF  Y*r  1 HE  N 19b 
I H b IF  Y * I 1HFN  19? 

I P 6 IF  Y*?  THEN  1 PH 

IPP  IFT  FC2 >*E C?  )\l FT  F C | > *F  C 2 > /T  2 ♦ F C 3 > 

I h 9 IFT  F C?  >«C  I -F  C?  > >/C  U *T  I MW4  »?*FC?> 

192  IM  F C 3 ) =V 4»2*FC? > ✓T?4W4t?*E( 3)-£ c I )/( lfr *T  I > 

191  GOTO  199 

19?  LET  F C2>=l -EC2I-FC? >\LET  F C I )■<  I-|C2)-E<2)>/12-EC I > -E < 3 > 

193  LEI  F (2  )*-EC  I )/TAM  ET  FC3>«2 

194  CJTO  199 

19b  l F 1 F C 2 > *F  f 2 > \l  FT  F C I ) *E C 2 ) /T 2* E C | > 

196  IM  FC?  >*V  4t?*C  I -EC?  > >4EC2>  XC  T2*  1 I >*EM  >/ 12 

197  IFT  F (3>*W4t?*( I ✓T2-EC3)-E<2)/T2)*EC I >✓< 12*1 I > 

199  DATA  "XC2)"#"X<  I >" # "X C 2 >" # "X C 3 >" 

21 b l ET  1*1 
???  F UT  P # 1 32  » I 
23b  LET  1*2 
? 42  F*L  1 H # b2  * I 
?b?  CHAIN  "|  I NFL  I " 

?62  IF  Y? * I 1HFN  2 7r 

?6b  FC2>*1-DM?>\ECI>*2N. C?)*DbC4) 

?66  F C3>*-AI  SCM?>WHI  47  S IF  Y *F  THEN  172 

267  Fh|Nl  "POUNPAHY  CONOITIONb  F Oh  ln-2  CONDUCTION  HF.it  I OP" 

?6P  Fh  INI  "71 *"l  7 I 

269  FhINT  "F ’ C 2 > *" | F C 2 > » "F C I >*"IEC  I > . #,E  C 2 > *M  I E < 2 > * " E C 3 > ••' I E C 3 > NGOl  0 172 
?72  F C2  )*|  ~Pb(A)NE  C I J *2  Nf  (2)*Db<  4) 

? 7?  FC3>*-APS«EC2) )/M(4)  \ IF  Y »2  THFN  172 

273  FhINI  ••FOUNOAhY  CONDITIONS  F Oh  Th-I  CONDUCTION  FFhlOP" 

274  Fh  INI  "7  » *,,I7  I 

27b  F’hINl  **E  C2  )•*•  IF  <2  )*,fF  C I )*•’ »FC  I ) ,”E  C ? )■••  IF  C?  ). "|  C 3 >■’•  IM  37NG010  172 
322  ENT 


•:i^Jtsrrrz:r.~  : ,r- 


1 


1 


? 
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PROGRAM  SEGMENT  "LINEQl 


OLD 

OLD  PHOCHAF  NANE--L  I NEl  1 
HEADY 


LIST 

tee 

i if' 

120 
1 22 
12S 
1 30 
I 3b 

1 40 
1 8 0 
1 90 
191 
1 92 
1 93 
1 94 
I 9 b 
1 96 
1 97 
1 9ft 

1 99 
20.0 
201 
PH? 

203 

204 

20b 

206 

2H7 

20ft 

209 

210 

21  I 
212 

213 

214 

2 I b 
216 
220 
221 
222 
223 

22  b 
226 
227 
22  ft 

230 

231 

232 

234 

236 

237 
23ft 
239 
2 50 
2b2 

2b4 

256 

257 
2A0 
290 
2 92 

293 

294 

29b 

296 

300 


Hfcft  PKOGhAF  10  SOLVE  4-L1NEAK  EQUATIONS 
HECOHD  AC3#3>#FC3>.,LC 4> 

KECOftD  X(3)#D4<20) 

OPEN  ft  # " DA T A 1 " 

LEI  1=1 
CET  ft# 1 10#  I 
1*2 

CEl  ft  # 1 20  # I 
PIP  P(23)#T(3#3) 

LEI  P<0)=AC2#2 )*A<  3#3)~A(2#3)*A( 3# 2 > 

LEI  P<  1 )=A<1#2)*AC3#3)-AC|#3)*A<3#2) 

LET  P (2 ) =A  < 1 #2)*AC2#3)-A(2#2)*A< I #3) 

LET  P<3)*A(2#1)4A(3#3)-A<2#3)*A«3#1 ) 

LF1  FC4)=A(1 #1 >*A<3#3>-A( 1,3)*A<3# I ) 

LET  P < b > = A ( I#1  )*AC2#3)-A(  ]#3)*A(2, 1 ) 

LEI  P(6)=ACP#1  )*A(3#2)-A(3#|)*A(2#2) 

LEI  ft  < 7 ) = A ( 1 #1  )*A(3,2)-A(3#1)*A(1#2) 

LET  ft<p  > = A( 1 # I >*AC2#2)-A(2# I )*A( 1 #2  ) 

LET  ft(9)=A(0#2)*A(3#3)-A(0#3)*A(3#2) 

LET  ft(^l0)=A(0#2)*A(2#3)-A<0#3)*A<2#2) 

LET  ft<l  I >=A<0, 1 )*A(3#3)~AfK,3)*Af  7*  1 ) 

LET  PC 12>*A<0,1  )*A<3,2)-A(0,2)*A<3,  I ) 

LET  F < 13)=AC0#2  ) • A ( 1 ,3)-A(0#3)*A<  1,2) 

LET  PCI 4)*A<0# I )*A(2,3)-AC0,3)*AC2, ) ) 

LET  PC  1 b)*A<0# 1 >*A(2#2)*A<0#2)*AC2# 1 ) 

LET  P ( 1 6>*A<0, 1 ) • A ( I ,3)-AC0,3)*A<  I# I ) 

LET  P< 1 7 >=AC0, I > • A ( | ,2)-ACP,2)*AC I#  1 ) 

LET  PCJft)=A<P#P>*A<  1,1  ) • A ( 1,0  )*AC0, 1 ) 

LET  ft(19)=AC2#0)*A(0,  1 )-A(0,0)*AC2,  1 ) 

LEI  P(?0)*A<0,0  )*A< 3, | ) -A ( 3# 0 )*A( 0, I ) 

LET  ft  C2  I )*ACI#0)*AC2# 1 )-AC2,0>*AC I, I ) 

LET  PC22)*A(3#0  )*A( 1,1  )-A(l,0)»AC3#l  ) 

LET  PC23)*AC?,0)*AC3,1  )-A(3,0)*A(2# I > 

PH  I NT 

LET  D*P  < 0 ) *P  ( 1ft  > ♦P  ( 1 >*PC  I 9>4ftC2>*ftC20)*ftC9)4pC2l  >4pC  I0)»ft<22) 
LET  D*D*P( 1 3>*PC23> 

LET  TC0#P)*CAC  I , 1 >*PC0)-AC2#  1 )*p(  1 ) ♦ A ( 3#  1 ) • ft  ( 2 ) )/D 
LET  TC 1 *0)  = C -AC 1 ,D)*H<0  ) ♦A(2,0)*P(  1 )-AC 3#0)*PC2)  )/D 
LET  TC?#0)=<A(l,0)*PC3)-A(2,0)*P(4)*A(3,0)*PCb))/O 
LET  TC3#0>=C-AC I #0>*R< 6>4AC2,0>*ftC 7)4 -AC 3,0>*ftC6>  >/D 
LET  T (0# 1 )*  < -A(0# 1 ) »ft  < 0 ) ♦AC?#  1 )*PC  9 >-AC  3#  I >*ft  C 10)  ) /D 
LFT  T C 1 , 1 ) = (A(0,0)*P(0)-A(2#0)*PC9)+A(3#0)4PC10))/D 
LET  T ( 2 # 1 ) = C-  A(0#0)*PC3)^A(?,0)»PC1I>-AC3#0)*PC14))/D 
LET  T ( 3 # 1 ) = (AC0,0  )*P<  6)-AC2#0)*ftC  12  ) ♦ A C 3#  0 ) • P < lb))/D 
LET  T(0,2)*(A(0#1  ) *p  C I )-A(  1,1  )*M9)*A(3»I  ) *ft  C 1 3 ) ) / D 
LET  T(1#2)  = (-AC0#0)*PCI  )4A<l#0)*ft<9)-A(3#0>4ftCI3))/D 
LET  TC2#2)*(A(0#0)*P(4)-A( 1 , 0 ) *ft  C I 1 >4AC3#0>*ftC16>)/D 
LET  T(3.2)*(-A(0,0)*p<7)4ACl,0)*ft(12)-A(3,0)*P(l7))/D 
LET  T(0,3)*C-A(0,l)*p<2)4A(l#l>*ftCI0)-A(2,|)*ft(13))/D 
LET  TCl#3)*(A(0#0)*p(2)-AC1,0)*PC10)4A<2#0)*b<13))/D 
LET  TC2*3)=C“AC0#0)*PCb)*A< I # 0 > • P C 14)-AC?,0)*PC 1 6 ) ) / D 
LET  T <3#3>*CAC0,0)*PC8)-AC |,0>*PC I b > *A C 2# 0 ) *P C 17) )/U 
LET  X<0)r7C0,0)*FC0)+TC0,l)4ECI)*TC0,2)*F<2)*T(0#3)*F(3) 

LET  Xfl  > = T ( l#0)*F(0)4T( 1# 1 > • F C 1 )♦!< 1,2)*F(2)*1< 1#3)*FC3) 

LET  X<2)*T(2,0)*FC0)4T(2,l)*F<l)4TC?,2)*F(2)4Tf2,3)*FC3) 

LET  XC3)  = TC3#P)4FC0>47f3#J  >*F(I)#TC3#2)*F(2)*TC3,3)*F(3) 

PHINT 

LET  1*2 

PUT  ft, 120* T 

IF  D4( 1 0 )*0  THEN  29b 

I F D4C 1 0 )* 1 THEN  296 

FW  I NT  ''SPECIFY  D4CIP)  FO  h CHAIN  S T A TEPE  N T'*  Nftft  1 NT  VG  01  0 300 
CHAIN  "EftftP 3" 

CHAIN  ••EKHPA** 

END 


PROGRAM  SEGMENT  "EKPP3" 


OLD 

OLD  PROCKAP  NAHE--EKPP2-3 
HEADY 


2b  OPEN  8 . "DA T A 1 ** 

30  RECORD  G<4).HC6>.Db(7>.K(4> 

32  KECOKD  X(3).D4<20> 

34  LET  I = 0 
36  GET  8.30. 1 
38  LET  1*2 

40  GET  8.32.1 

bP  PKINT  "04(2  )*••  JD4<2>,"L*<  6>*"  JP4(  6) 

1 34  I b*K(0 >\7b*K( I >\T6*K(2>\72*KC3> 

lb?  Y*Db(0>\PR  INI  ,,Y*MI  YNPKIN1  ” Tb«" J T bNPK l NT 

154  7 1 *Db( b >NPK I NT  "7 1 *" I 7 1 

I b9  73*71  -4*INT( 71 /4>nYI «C - I ) *INT( 73)NY2*C - I ) «INT<  Z3/2) 

160  IF  Y I * I THEN  162 

161  IF  Y 1 * • 1 THEN  164 
16?  IF  Y?  * I THEN  166 

163  IF  Y2  * * 1 THEN  167 

164  IF  Y2  * - 1 THEN  169 
16b  IF  Y2  * 1 THEN  168 

166  PRINT  "TKANSIS10K  1 CONDUCTION  PEK I OD*' NCOTO  171 

167  PKINT  "TRANSISTOR  2 CONDUCTION  PE  K I OD”  NGOTO  171 

168  PKINT  "DIODE  I CONDUCTION  PEKIOD"\GOTO  170 

169  PKINT  "DIODE  2 CONDUCTION  PEKIOD"NGOTO  170 

170  GOTO  172 

171  72*1 

17?  8 * SI  K f ( H ( I >♦*<  3 >/X(2>  ) *2*H<  2 ) *2  >*SGN(HC  1 > *X < 3 > / X ( 2 > > 

173  Nl«7?*3.1 41 b93/M<2>\N2«N| /2SV0IP  V0(3>.J0C1> 

174  IF  Y * 0 THFN  I 7 bNGOT  0 438 

I 7b  D0*APSC2*N?*CC0 )4HC2 ) »2*XC? ) )NPK INT  "D0*" I D0NP* I NT  "NI*"INI 

178  PKINT  "PlTP”NPK|NT 

180  FOK  12*0  TO  Nl  SI EP  N2 

18b  I FT  7*H<?)*T?*ATN(HC?)/(H( I )4XC3)/X<?)>) 

190  l El  V0(1)*X(0)4XC1  )*tXPCH(0)*T2)*XC2)*EXPCHC  I > • T2  > •*  * S 1 N ( 7 > /H  < 2 > 
19b  LET  J0C  I )*C(0)*(X(  | )*H( 0 )*EXPCH( 0)*T2 ) ) 

196  LET  J0C I >*J0( | )4G<0)*<H< I )*XC2)4EXP(H( I )*T2)*8*SIN<7)/HC2>) 

107  LET  JPM  >*J0M  >*G<0 >»<H<2>*X(2>*EXP<H<  1 >» T2 >*K*COS< 7 >/H(2 > > 

198  IF  T2  «•  34N2  THEN  260 

199  IF  APSfJ0M>>  «*  D0  THEN  20b 

200  GOTO  260 

21/  b LET  NI*T2\PKINT  "Nl  *"  IN  1 \PK  I NT 
?|0  IF  Y I * - 1 THEN  220 

211  IF  J0< I )«0  THEN  213 

212  D4Cb>*V0C 1 )ND4C6)*J0C 1 JNCOTO  ?37 

213  J0C  I )*0W0(  1 )*D4(  b)\D4C  6)*J0C  I ) 

21b  GOTO  ?3 7 

220  IF  J0<l>»0  THEN  2?5 

???  D4< 4>*T?ND4<0 >*V0C I )ND4(2)«J0C | >\COTO  237 
22 b J0C | )*0\V0C I )*D4(0)\D4(2)sJ0( | ) 

237  IF  Y2  * 1 THEN  ? 40 

239  V/0  ( I ) * I - V0C  I >\J0(1  >*-J0C|  ) 

? 40  PKINT  INTC9999*(TS*T2>  'T6>. I NT < S00P* 4999* 90 < I >/2.b> 

24?  PKINT  |NTC9999*CTb4T2)/T6)."  "I 

?4b  PKINT  I NT  C b000 ♦ 4999* J0 ( I > /0  • | > \ PH  I NT 

Pb0  GOTO  300 

260  IF  Y I «-l  THEN  26b 

P 6 1 IF  JP<l>«0  THFN  ?68 

?64  D4C  b)*V0(  I )ND4C  6)*J0C  I >\GOTO  •>  70 

?6b  IF  J0 < 1 ' *0  THEN  267 

? 66  PM  4 ) r T2\D4(0  )«  V0(  I )M  MP)*J0(  | )\  rCTO  270 

267  J0(  1 )*0Wf.C  I )*D4(0)\D4C?)*J0f  I '\rr»To  270 

268  J0C  I )*0\V0<  1 )«D4C  b)ND4C6)*J»f  MNGOTO  270 
?70  IF  Y?  * I THEN  280 

27b  I >*l -V0C I )NJ0C 1 )*-J0C | ) 

280  PRINT  INI C9999*(Tb«T2>/T6). !N1Cb000«4999*V0( I )/2.b) 

2Pb  PMVT  !N?<9999*CTS4T2>/T6>."  "I 

P86  PRINT  INT(b000*4999*J0< I > /0  • | >\ PRINT 
290  NFXT  T? 

.700  PK1NTNFKINT  "PI  TT"NPR  I NT 


« 


I 


I 


I 

* 


4l 


3?b  L FT  D?=D P/2P 

3P6  PRINT  "D4C2>="JD4<2>,"D4<6)="lD4< 6> 

3?7  IF  72* 1 THEN  396 

3?R  PRINTSPnINT  "KLTP"SRhINl 

3ltf  FOR  T2*N  1 TO  1.1P*NI  ST  FP  N'2/2? 

31b  LET  7sHC2)*T?4ATNCMC2)/CHC  1 )♦>! 3)/>C2>)> 

32?  LET  VP<I)=X<?)  <( I >*EXR<H(?)*T2>*XC2)*EXP<H< I >*T2)*K*S1NC? >/H<2> 

33?  LFT  JPO  )=C(P>*  <X<  I > *H  < fl  > * E X P < H < 0 > * T2  > ) 

34?  LET  JPf  I >= JP<  I ) *C(  V > • ( H < 1 >*XC2>*FXP(MM  >*T2>*X*SIM7)/M<?>> 

3 S?  LET  JPC  1 >=JPC  I >4C<?)*<H<?>*X<?)*EXP<H<I  >*T2  >*K*C0Sf 7 >/M<2 > > 

362  IF  APS<J0C1>>  « = PP  THE  N 364 

363  GOTO  3P? 

364  IF  Y I =•  I THEN  36f 

365  IF  D4C6)*?  THEN  3 7PSP  4(  b ) = VC  < I > \ C OT  0 372 

366  IF  D 4<2>=P  THEN  36P 

367  D4( 4>=T?SD4<? )=VP< I ISCOTO  372 

36P  JPf  1 )*PWf(  1 >=D4(W)\U4C?>*J?(  | ISGOTO  372 
37?  JPf  1 >=?SV?f  I )=D4f  b)Sp4f  6 >*JPf  I ) 

372  IF  Y?  = I THEN  376 

374  VPf  I ) = I - VPf  I ) SJPf  I >*-JPf  I ) 

376  PRINT  INT  f 9999*f  Tb-*T2  >/76>,  iNlf  SP0?*4999*VPf  1 >/?.b> 

377  PkINT  INTC9999*CTb4T2) /T6)#,f  "I 

378  PkINT  INTfSPP?44999*JPf 1 )✓?. I >SPR |Nl 

379  LET  Nl  *T2\FK  INTNFR  INT  "N  I «"  l N I \PR  I M SCO  T 0 396 
38?  IF  Yls-I  THEN  383 

3M  IP  C4<6>=0  THEN  38b 
IP?  P4(  5)  = V?f  I >\GOTO  3P6 

3P3  IP  P4f2>=?  THEN  3H 4\p 4C 4 ) = 1 2 M>4f ? ) * VP f I > SGO T 0 386 
3P4  JPf  I )*?SVP<  I >=l)4(p)\r)4(2>sJP<  I I^GOTO  386 
3PS  JPf I >*?SVPf I >=D4f b)SP4f 6>=JPf  I ) 

3P6  IP  Y? = 1 THEN  39? 

3PP  VPf  I )s| -VPf  I >SJPf  ! >*-J?f  I ) 

39?  PR I N T INTC9999*(Tb*T2)/T6)#  1NTC 5P??44999*V?C I )/2.b) 

39?  PRINT  INT  ( 9999 • f T b* T? > /T 6 > # " "l 

394  PRINT  INT(  b?P('*4999*JPf  I ) /P  . I )NHb  INT 
39b  NEXT  7? 

396  72  *8  f 3 > 

39P  PRINT  SPRINT  "PL  I TMSPrt  INTS  I F 7?*P  THEN  4?b 
399  IF  Y2  * I THEN  4? 5 
4?P  IF  Y?  * • I THEN  422 
4?b  IF  Y 1 * 1 THEN  42? 

4 1 5 Dbf6)*VPfl  )NDbC  7 > * jp  f | ) spbf I ) =N  I \ GOTO  429 
42?  Dbf2)*V?f I >SPbf3)=J?f I > Spbf I > *N  I \ GOTO  429 

422  IF  Y I * I THEN  42  4 

423  Dbf  6 > *1 -90(1  )\DS<  7 ) *- JPf | ) \ pbf I )»N1  S GOTO  429 

424  Db(2)*l -VPf 1 >spb(3)=-J0f 1 ISDbf 1 >*NI  \ GOTO  429 

425  IF  Y2  * - 1 THEN  42  7 

426  Pbf2 >*VPf I >\P5< I )=NI SGOTO  429 

427  Dbf 6>*l -VPf 1 >SDbf 1 > =NI \COTO  429 

429  Y * Y ♦ I \ D5C?)*Y  \ PRINT  "Pbf  I > s" I Dbf  I )\PR|NT 
43?  LFT  1 =P 

432  PL 1 T 8#3P# I 

433  1=2 

434  PUT  8,32*  I 

435  GOTO  49b 
43H  72*K<3> 

439  PRINT  " D4  f 2 )*"|P4(? )»"P4( 6) *" J D4f 6) 

44?  LET  NI*Pb(l  ISLET  N?  *N  1 /25SPR I NT  **NI  *"  IN  I SpR|NTS|F  Y*1  THEN  472 
44b  IF  Y * 2 THEN  4b2 

446  PRINT  "ERROR  IN  ASICNPENT  OF  Y " S GOTO  6?? 

4b?  PR  1 NT SpR 1 NT  "PI  TP"SPR1NT 

4b3  PRINTSPRINT  "7  PE GRFF" I " TRANSF.  VOlTAtF"#"  OUTPUT  VOLTAGE" 

4bb  FOR  T? * ? TO  N I • I . PPb  STEP  N2 

4b6  LFT  7*Mf2>*T?4ATNfMf?)/fMf  I )*XC3)/XC2>)> 

4b7  LET  V? f2)=X(P)#Xf  I >*FXP(M<P>*T2>*Xf2>*|XKH<  I >«T2)*K*S|N<?>/Hf2> 
4b#  LFT  VPf 3>=APSf V0f2>>/Gf 4>SIF  Y1*-l  THEN  46? 

4b9  J?C 1 )=D4C6)SIF  J0f|>*P  THEN  46?  SD  4t  8 ) * VP  f 2 > SG  0 T 0 463 
46?  JPf  I ) sD4(? )S | F JP  Cl)*?  TMFN  46? 

461  P4C3)*V?(?)SG0T0  463 

46?  V?C2 ) *-D5( 4 )*FXK - ( T2 -p4( 4) ) /Hf  4) )\ V?( 3 ) «APS( V?<2 ) ) /G<  4) 

463  IF  Y?  * I TMFN  46b 

4 64  V(  (?>=-VPf2  > 

46b  PKINT  INTfP999*(Tb4T?)/T6)»  I N T f bPPP ♦ 4999* VP f 2 > / I •?) 

466  PRINT  INTC99994Y Tb4T?)/T6)*"  " » I NT C bPPP* 4999* VPf 3 ) ✓ b > SPK | NT 

46P  NFXT  T?S  GOTO  484 

47?  PRINTSPRINT  "PVTL"SPRINT 


473  HUNT  "7* DEGREES"," INDUCTOR  VOL  TACE"\EOK  T2*0  TO  NI*l.00S  STEP  N2 

474  LET  7=H(2)»T2+ATN<Ht2)/<H<  I ' + x< 3 > /X « 2 > ) > 

47  5 LET  V0<0>*X<0>+X<l>*EXP<H<0>*T2>+X<2)»EXP<M<n*T2>*K»SlN<?)/H<2> 

476  IF  Y 1 = - I THFN  478 

477  J0(  I ) =D4<  6 > ' I F J0<l>=0  THEN  5 1 PND4<  7 ) * V0<  0 > \GOTO  460 

478  J0< I ) *04(2 ) VIE  JP<l>=0  THEN  479\D4( I ) * V0 1 0 ) \GOTO  480 

479  V0<P)=P\D4< I )=0 

480  IE  Y?  = I THEN  482 

481  VP(P>=-V0<0> 

48?  PRINT  INT<9999*<T5+T2)/T6>» I NT< 5000+ 4999* V0 ( 0 > /2 . 5 > VPk I N1 

483  NEXT  T2\Y=Y+I \D5< 0 > =Y\GOT 0 490 

484  IE  Y2  = I THEN  486 

485  V0<2)=-V0<2> 

486  P5< 4) = V0<? )\Y=0\D5(0)=Y  \ IE  72*0  THEN  488 

487  7 I =7  I ♦ I \GOTO  489 

488  71=71+2 

489  D5<  5>=7I \T5=T5+D5<  I )\M0»=T5MF  7I»7S  THEN  496 

490  PRINTNPRINT  "PL  TT"\ I = 0\PUT  8»30>I 
49?  LET  I =2 

493  PUT  8.32. I 

494  IE  7 1 +2*75  THEN  520 

495  CHAIN  "Exm" 

496  1*2 

497  PUT  8,32,1 

498  I=0\PUT  8.30,  I 

499  GOTO  600 

510  V0<0)=0\D4< 7)=0\GOTO  480 
520  D4« 10)=0\GOTO  495 
600  PRINT  "PLTT"VEND 


^3 


PROGRAM  SEGMENT  "EKD1 


OLD 

OLD  PROGRAM  NAME  — FK D 1 


HEADY 


LIST 

50  PRINT  "PROGRAM  TO  PRINT  CONTENT  OF  DATA  1 ! " 

60  PRINT 
70  PRINT 

100  RECORD  G< 4) ,H<6),D5< 7>,K< 4) 

200  RECORD  A(3,3),FC3),L<4) 

300  RECORD  X(3),D4(20> 

400  OPEN  8 , "DATA  1 " 

420  LET  1=0 

430  GET  8,100,1 

440  LET  1*1 

450  GET  8,200,  I 

460  LET  1=2 

470  GET  8,300,  I 

500  FOR  1=0  TO  4 

510  PR  INT  ”G<  " » I J ” ) =" » G(  I > , 

512  NEXT  I \IRINT 

515  PRINT 

520  FOR  1=0  TO  6 

530  PRINT  "H("J I }">="JH<  1 ) , 

532  NEXT  I SPRINT 

535  PRINT 

540  FOR  1=0  TO  7 

550  PRINT  "D5C"! I J" >=";D5( I ), 

551  NEXT  I SPRINTSPRINT 

552  FOR  1=0  TO  4 

553  PRINT  ”K  ("  ) I J " ) =" I K ( I ), 

554  NEXT  I SPRINTSPRINT 

555  PRINT  "MATRIX  A<3,3>" 

560  FOR  1=0  TO  3 

565  FOR  J=0  TO  3 

566  PRINT  A ( I, J), 

567  NEXT  JSPRINT  " " 

568  NEXT  I SPRINT 

569  PRINT 

570  FOR  1=0  TO  3 

575  PRINT  "F<" J I J" >=" JF< I >, 

576  NEXT  I SPRINTSPRINT 

577  FOR  1=0  TO  4 

578  PRINT  "L("J I J")="JL<  I >, 

579  NEXT  I SPRINTSPRINT 

580  FOR  1=0  TO  3 

58  5 PRINT  "X ( " J 1 I " ) =" I X ( I ) , 

586  NEXT  I SPRINT 

588  PRINT 

590  FOR  1=0  TO  20 

595  PRINT  "D4("J I J")="JD4( I ), 

596  NEXT  I SPRINT 

597  IF  04(10)=!  THEN  599 

598  COTO  600 

599  CHAIN  "EKD2" 

600  END 
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PROGRAM  SD3MENT  "EKD2" 


1 


OLD 

OLD  PROGRAM  NAME--EKD2 

READY 

LIST 

40  PRINT  "DATA  RECORD  CHANGE” 

50  RECORD  G(4)#H(6)#D5(7)#K<4) 

55  RECORD  X<3)#D4<20> 

60  OPEN  8# "DATA  1” 

61  FOR  1=0  TO  4YINPUT  G<I)\NEXT  I 

62  FOR  1=0  TO  6\ INPUT  H(I)\NEXT  I 

63  FOR  1=0  TO  7 \ INPUT  D5(I)\NEXT  I 

64  FOR  1=0  TO  4 \ INPUT  KCI)  \ NEXT  1 

65  FOR  J=0  TO  3 \ INPUT  X(I)  \ NEXT  I 

66  FOR  1=0  TO  6 \ INPUT  D4( 1 > \ NEXT  I 

67  FvOR  1=7  TO  13  \ INPUT  D4(I)  \ NEXT  I 

68  FOR  1=  14  TO  20  \ INPUT  D4( I ) \ NEXT  I 
70  1=0  \ PUT  8/50#  I 

80  1=2  \PUT  8# 55#  I 

85  IF  CD5(5)+20)=K( 1 ) THEN  95 

90  GOTO  100 

95  CHAIN  "EKPP1” 

100  END 
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